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induced by recurrence properties of random walks 
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Abstract: We study a quantum walk (QW) whose time evolution is induced by a random walk (RW) 
first introduced by Szegedy (2004). We focus on a relation between recurrent properties of the RW and 
localization of the corresponding QW. We find the following two fundamental derivations of localization of 
the QW. The first one is the set of all the I 2 summable eigenvectors of the corresponding RW. The second one 
is the orthogonal complement, whose eigenvalues are ±1, of the subspace induced by the RW. In particular, 
as a consequence, for an infinite half line, we show that localization of the QW can be ensured by the positive 
recurrence of the corresponding RWs, and also that the existence of only one self loop affects localization 
properties. 

CT 1 1 Introduction 



A frame work for quantum search algorithms by using the coin-shift QW on a spatial structure was con- 
structed by [TJ [21 13] for example. Its one step time evolution U is given by U = SC: the first operator C 
CN| . is so called coin flip assigning a unitary operator called quantum coin (in many case, the Grover operator) 
to each vertex, where the dimension of the quantum coin is the degree of the corresponding vertex, and the 
second operator S is called shift which moves a quantum particle to its neighbor. The paper 4 is one of 
the nice reviews on these works. Szegedy [5] proposed a generalized QW whose time evolution comes from 
a corresponding RW on a bipartite graph. This walk is called bipartite walk or the Szegedy walk. Let the 
time evolution of the Szegedy walk and a coin-shift QW be W, and U, respectively. Then as mentioned by 
[6] for example, an appropriate choice of the coin flip operator of U gives W — U 2 . The eigensystem of W 
is explicitly obtained in [5] ■ In [El H] , they blush up the quantum search algorithm and a notion of quantum 
hitting time of [5] with the walk. 
^ ■ In these works, the eigensystem of the subspace induced by the spectrum of the corresponding RW, 

named H^ R ' in this paper, plays an important role: the spectral gap of the corresponding RW determines 
the quality of the quantum search algorithm. Moreover graphs treated in these studies are finite because 
of its algorithmic application. Since is invariant under the action of U and W , an appropriate initial 

state gives us considerations outside of its orthogonal complement space %^ R > = ri^ s \ In this paper, we 
take an initial state which is a convex combination of u g rl^ and v € 'HS S > to get some effects of the 
subspace rl^ on behaviors of the QW, and also treat a simple infinite graph, that is, a half line with self 
loops. We focus on the time averaged limit measure. Due to its infiniteness of the system, the distributivity 
of the time averaged limit measure is not ensured. In fact, except in the case of positive recurrent case on 
the half line, the time averaged limit measure of the RW becomes the null measure. In this paper, if we 
can take an appropriate initial state so that the time averaged limit measure of the QW is positive, then 
we say that localization of the QW occurs. We show that all the £ 2 summable eigenvectors of the RW and 
the orthogonal projection onto rl^ s ^ of the initial state are essential to the derivations of localization of the 
QW. Off course, the existence of a finite number of self loops has little influence to the behavior of the RW. 
However as we will see, since the number of the self loops determines the dimension of the subspace of T-L^ 
in the QW, the behavior, especially localization property, is quite sensitive to the self loops. Actually, adding 
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only one self loop to even a transient RW gives a change from a QW without localization to a QW which 
exhibits localization. 

This paper is organized as follows. In Sect. 2, we give a review on the Szegedy walk and clarify the 
eigensystem of the time evolution U on graph G(V,E), where V is the set of vertices, and E is the set of 
edges. This paper treats graphs with self loops. We denote the set of of self loops as5C-E. Since the time 
evolution of [5] W, is expressed by W — U 2 , the statement of part 4 in Theorem 1 of [5] implies that the 
possibilities of eigenvalues of are {1}, {—1}, and {1,-1}. In this section, we find that the eigenvalues 
of are ±1 with multiplicities \E\ - \V\ + m(l) - \S\ and \E\ ~ \V\ + m(-l), respectively. Here m(±l) 
are the multiplicities of the eigenvalues ±1 of the corresponding RW. Section 3 proposes a relation between 
recurrence property of the RW and localization property of the QW. It is obtained that the time averaged 
limit measure of the QW is bounded below by the stationary distribution of the positive recurrent RW. We 
also see that this localization comes from HS R \ On the other hand, we find that the effect of TiS^ to the 
localization property on an infinite half line is determined by the number of self loops and the recurrence 
property of the corresponding RW. More precisely, the support of localization in {0, 1,2,...} with a special 
initial state depends on the recurrence property of the corresponding RW with finite self loops. Finally, we 
give a discussion. 

2 Review on the Szegedy walk 
2.1 Definition of the quantum walk 

In this paper, we denote the Hilbert space of square-summable functions on a countable set X by £ 2 (X). 
The inner product is defined by 



and the norm by ||/|| = (/, f) 1 ^ 2 - We define a complete orthonormal system of £ 2 (X) by {S x } xe x with 



The graph G(V, E) treated here is connected and all the degrees of vertices are finite. In this paper, we treat 
graphs with self loops. The self loop means (u, v) G A(G') with u — v. We denote S as the set of all the 
self loops. At first, we consider an ergodic RW on the graph G(V, E). Let M be the stochastic matrix of the 
RW, (M) with J2uee{G) ~ 1- ( a wai ker located in position v G V(G) jumps to the neighbor u 

with probability p u>v . ) We also give a self adjoint matrix J defined by (J) u ,v — y/Pu,vPv.u which plays an 
important role in this paper. Recall that if M has a reversible measure 7r(j) = ttj, then PijiTj — PjiiTi- This 
relation is called "detailed balance condition". When ttj ^ for every j, then J = D-^MD 1 / 2 , where 
I? is a diagonal matrix with (D)jj = ttj. In the following, we should review recurrence properties of RW, 
which appear again to the localization of QW. For more details, see [HI [TU] for example. If the probability 
that a number of revisits to the starting place is infinite equals to 1, then we say that the walk is recurrent, 
otherwise transient. Moreover if the mean first hitting time in the recurrent walk is finite, then we call it 
positive recurrent, otherwise null recurrent. 

The symmetric oriented graph of G(V, E) is a digraph G'(V, A) whose vertices and arcs are defined as follows: 
V(G') = V(G) and {it, v} G E(G) iff (u, v), (v, u) e A(G'), where (u, v) is the arc from the vertices v G V(G') 
to u G V(G'), From now on, we give a QW corresponding to the above RW. The total space of the QW is 
H = £ 2 (A) whose orthonormal basis are associated with the arcs of its symmetric oriented graph G', that is, 



Note that the dimension of our quantum walk is \A\ = 2\E\ — \S\. The time evolution U — SC is defined as 
follows: 



{f,g) = ^2f(x)g(x), f,ge£ 2 (x) 




n = 1 2 {A) = span{ : («, «) G A{G')}. 
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• Shift S: For any £(„,„) with (u,v) e A(G'), 

<S , «5( T4 ,«) = (1) 
This is a permutation on A(G') called "flip flop" [2:. 

• Coin flip C: Let a u and b u (u £ V(G')) be 

Q"U — ^ ^ -\/Pv~u^ (^v .u) s bit — SCL U . 

v:(v,u)eA(G') 

We call a u and b u "incidence vector of it" and "swapped incidence vector of it", respectively in this 
paper. Denote Ha and IIb as the orthogonal projections onto the subspaces of Ha — span{a M : u G 
V(G')} and Hb = span{6 u : u € V(G)}, respectively. Then the coin flip C is defined by 

C = 2n A - /. (2) 

Noting S 2 = /, the one-step Szegedy walk, W, is equivalent toW = U 2 — ref^refA, where ref^ = 2Ha — I(= 
C) and refs = 211b — /. An equivalent expression for the QW is that: a quantum particle on v € V at 
the present time step which was in v € V at the previous time step, moves to its neighbor, w € V , with 
amplitude 1^/p w ,vPu,v — 3w,u at the next time step. So only if a quantum particle returns back to the same 
place, an extra term — 1 appears in the amplitude associated with its moving. When the corresponding RW 
is symmetric in that p UyV = l/deg(w) for every v € V, where deg(w) is the degree of v, the QW becomes so 
called the Grover walk which has been intensively investigated. 

For \&o € £ 2 (A) with H^oll 2 = 1) let b e a random variable which is determined as follows. 

p(x^ o) = u) = K*(»,u), ^**o>| a . 

v:(v,u)eA(G) 

The random variable A^* -* corresponds to a position when the QW with the initial state V^o is measured 
at time t. 

2.2 Eigensystem of the QW 

Let A be a weighted incidence matrix of G defined by 

u£V{G) 

Also define B as 

B= 

uGV(G) 

Remark 1 The eigenvalues of the shift operator S are ±1 with multiplicities \E\, and \E\ — \S\, respectively. 
The eignespaces of eigenvalues 1 and — f for S are described by 

M+ = span {{ 8 {u ^ v) + : (u, v) E A(G') with u^v}U {<$(„,„) : (u, v) E S}} , (3) 

M- = sp&n{5 {UjV) - 5( VtV ) : (it, v) € A(G') with u ^ v}, (4) 

respectively. 

The "correspondence" between the RW and the QW derives from the following lemma. 
Lemma 1 
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1. For x e [—1,1], let p{x) G P(V) satisfy the following eigenequation: 

Jp(x) = xp(x). (5) 

Then we have 

Uq ± (x) = e ± ^ q± (x), (6) 
where x = cos9(x), and q±(x) e £ 2 (A) with 

j(I-e±^S)Ap(x) :x${±l}, 
q±[) ~\Ap(x) :xe{±l}. [) 

2. If the orthogonal complement of = span{'M.A,'HB} exists, then rl^ = iJ-L^) 1 - is orthogonally 
decomposed into invariant subspaces with respect to the action of U , such that = W.^ © 
where for all s± € 

Us + = s + , Us- = — s_ 
with multiplicities \E\ — \V\ + m(— 1) — |<S| and \E\ — \V\ + m(l), respectively. 

Remark 2 If \V\ < oo, then span{a u : u e V(G)} = span{Ap(A) : A € spec(J)} since (a u ,a v ) = S UiV , 
(b u ,b v ) — S UiV and from its finiteness of the system size, the eigenvectors of J give the complete system of 
£ 2 (V). Therefore we have 

n (R) = span{q ± (A) : A e spcc(J)}. 

All the real parts of the eigenvalues of VS R ^ are all the eigenvalues of J. On the other hand, the rest of 
eigenvalues are given by 1 and —1 with the multiplicities \E\ — \V\ — \S\ + m(l) and \E\ — \V\ + m(— 1), 
respectively. 

In the graph without both cycles nor self loops, then the dimension of our quantum walk is reduced to 
d = 2\V\ — 2. It is known that whe the graph is a path, then the spectrum of J, spec(J), has eigenvalue 1 
without multiplicity. If — 1 ^ spec (J), then from the above lemma, d = 2\V\ — 1 which contradicts to the 
dimension of the walk. Thus -1 g spcc(J). 

Remark 3 If the graph is a path, then TL^ =0. 

Let P{x) satisfy the following eigenequation. 

MP(x) = xP(x). 

The reversible Markov chain implies p(x) = D~ x / 2 P{x) which provides following remark: 

Remark 4 If the Markov chain is reversible whose limit measure ir(u) > for all u € V(G), then we can 
reexpress the eigenvectors of U as 



q±(z) = 



'd-V2(/ - e ±ie ^S)AP(x) : x £ {±1}, 
D-^ 2 AP(x) : x e {±1}. 



Now, we give the proof of Theorem 1. 

Proof. It should be noticed that AA^ = U A , BB^ = U B . Then wc have 

Ap(x) 4 Ap(x) A Bp(x). 
On the other hand, from A'B = B^A = J and the definition of p(x), 

Bp(x) m> 2xAp(x) - Bp{x) A 2xBp(x) - Ap(x). 
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Therefore the subspace spanned by Ap{x) and Bp(x) is invariant under the action of U. Thus we arrive at 

the conclusion of part (1) of the lemma. 

It is noticed that s G implies Ss G "H (s ' ) . We find that 

s h> -s H- — Ss, Ss i-4 —Ss A- — s. 
Therefore U(I ± S)s = =f(I ± S)s, for all s G H (s) . Since for any s G ((/ + S)s + (I - S)a)/2 = s, 

H+ = span{(7 + S)s; s G H (s) }, H- = span{(7 - S)s; s G H (s) }. 



Note that 

S(J + S)s = (I + S)s, S(I — S)s = —(I — S)s. 
Equations J5]) and <|9j) implies that 

n + = n (s) nM_, H- =u {s) dm+. 



Since 



S (j4j>(a:) + Bp(x)) = (Ap(x) ± Bp (as)) , 

5 (Aj>(a:) - BjJ(») = - (Ap(x) - Bp(x)) , for x G spec( J) \ {±1}, 
SAp(l) = Ap(l), for 1 G spec(J), 
SAp{-l) = -Ap(-1), for - 1 G spec(J), 



it is hold that 



dim (u [R) n M+^j = \V\ - m(-l), dim (u (R) n M-) = | V| - m(l). 
Combining Remark [T] with it, we obtain 

dim(K^ S) ) = dim(-H (s) n M-) = dim(X_) - dim{H {R) n M-) 



= \E\-\S\-{\V\-™{1)). 



In the same way, 



dim{H {S) ) = \E\-(\V\-m(-l)). 
Then we get the desired conclusion of part (2) of the lemma. 



(8) 
(9) 

(10) 



□ 



3 Sufficient conditions of localizations given by recurrence prop- 
erties of RWs 

3.1 Time averaged limit measure of QW 

The set of time averaged limit measures Moo is defined as follows: for any u,v G V(G), 



T-l 



Modu, v) = \ lim I £ P(X^ o) = u) : *o G J £ a w S {w , v) : £ 



\a„, r = 1 



This measure is so called Cesaro summation of the finding probability from v with the initial state \I/o to u. 
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Lemma 2 For any (u,v) G V X V, the time averaged limit measure ^^°\u,v) G M.oo(u,v) with the initial 
state = X)uj (u) v) a w$(w,v) satisfies the following statements. If the corresponding Markov chain has the 
reversible distribution ir, then the time averaged limit measure Jl^"\u,v) is bounded below as follows: 

-p£°\u 7 v)>\(a v ,* )\ 2 n(u)7r(v)+ £ |(<5 Ku) ,n 5 * )| 2 , (11) 
where Us is the orthogonal projection onto 

Proof. We consider the QW on G(V, E) with \ V\ — oo. The time evolution U is induced by the Jacobi matrix 
J. Let p{x) and q±(x) satisfy the eigenequations Jp(x) — xp(x) and Uq±(x) = a;q ± (:r), respectively. Fix 
a vertex vq € V called the origin and let us regard G as the graph G n — G n (V n , E n ) in the limit of n — > oo, 
where V n = {v £ V : p(vo, v) < n}, and E n = {{u, v} e E : u, v e V n }. Here p(u, v) is the shortest distance 
between u and u. Put the time evolution of the QW on G n as U n , and J n as the corresponding Jacobi 
matrix. Assume that the time evolution U n is decomposed into 

where is the normalized eigenvector of eigenvalue e 10 ^ . Lemma[T]and ||q±^(A)|| 2 = 2(1 — A 2 ) (A G 
spec(J„) \ {±1}), = 1 (A G spec(J„) n {±1}) provide 

{v^h D | (^== + %M>) ^ W ■■ A € spec(J„)| . 

We should remark that lim T ->oo Y^=q e it(e < jT = 5 Um . The time averaged limit measure fjS°' n \u, v) 

with the initial state \&o = J2 w -(w v)eA a w8( w ,v) is reexpressed by 

A& 0,n) K«) = £ £ l<«L B Uo>l>£ B \Wl a - (12) 

fc m:(tu,u)GA„ 

In the limit of n — > oo, put 

B = {x G [-1, 1] : p(x ) G £ 2 (V) \ {0}} (13) 

as the set of the mass points of J. Equation ([7]) in Lemma [T] implies that if P(xq) 6 £ 2 (V), then q±(xo) G 
^ 2 (A). Let 7ioo°' l (w, w) = lim n ^ 00 /lio' '"'(M, u). Then we have 

?i£ o) («,«)= £ £ ( VZ-^y + %o|=i>) |(q.(^o)^o)| 2 |(q £ (a:o)^ (tu , Il) )| 2 

+ E l(<W),n 5 * )| 2 . (14) 

w:(w,u)£A 

Since the corresponding Markov chain has a reversible distribution, combining Lemma Q] and Remark [2] with 
Eq. (JT3)l implies 

-p£°\u,v)> \(a v ^ )\ 2 7v(u)7v(v)+ E l(<W),n 5 * )| 2 . (15) 

□ 

Remark 5 J/|V| < oo, £/ien the time averaged limit measure is a distribution in that^2 uev ~jl^ \u,v) = 1. 

On i/ie other hand, if \V\ — oo, then its distributivity is not ensured in that X^«ev ~fi°°°\ u i v) < 1 as is seen 
in J2J/, for example. 
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Transient 


Null recurrent 


Positive recurrent 


Recurrence probability 


< 1 


= 1 


Ct 


< oo 


= oo 


Mean recurrence time 


= oo 


< oo 


C B . 


= oo 


< oo 



Table 1: Recurrence properties of RW 



3.2 Localization of QW on half line 

In this subsection, we treat a QW on a simple graph, semi infinite line with some self loops. We say 
localization happens at position i 6 {0, 1,2,...} when the time averaged limit measure with the initial state 
*o = Y,k-.\k-j\<i a kjti(k,j) satisfies 

Vg o) (i,j)>0. 

The stochastic matrix of the RW which induces the corresponding QW is given by 

\ 



M = 



f r ° 


<h 





Pa 


n 


12 





pi 









P2 





(13 







V 



where qj + rj + pj = 1, pj > (j > 0), qj > (j > 1). A particle at position j jumps to left, right and itself 
with probabilities qj, pj and rj, respectively at each time, independently. Thus the symmetric oriented graph 
G'(V,A) treated in this subsection is V(G') = {0,1,2,...}, A(G') = : \i-j\ = l}U{(i,i) : r, ^ 0}. At 

first, we review the recurrence properties of the RW on a half line. For more details, see [10], for example. 
Define 



Cn 



\ - qi ■ 



•9i 



Ct 



E- 



qi ■ ■ ■ qj 



If Ct = oo, then the RW is transient, otherwise recurrent. Moreover if Cr = oo in the recurrent RW, then 
the RW is null recurrent, otherwise positive recurrent. See also Table [T] We see the following interesting 
example given by [TU]: In the space homogeneous case pj = p, and qj — q (j > 1), then the RW is transient 
if p > q, if p = q, null recurrent, and if p < q, positive recurrent. On the other hand, in the following three 
cases (a), (b) and (c), we see that p^ = q^ in the limit ofi—¥ oo. However, we find that the three cases have 
different recurrence properties: case (a) is transient, while cases (b) and (c) are null and positive recurrent, 
respectively. 

(a) Pl = (i + 2)/(2i + 2), q z = i/(2i + 2), 

(b) Pl = (i + l)/(2t + 1), q t = i/(2i + 1), 

(c) p Q =l tPl = qi = 1/2, ^ = (i - l)/2i, q % = (i + l)/2i, (i > 2). 

Even if we add some finite number of self loops S = { j\ , fa , . . . , j n } C Z + to the above RWs so that pj + qj < 1 
for every j € S, there is no influence to these recurrence properties. This is a trivial statement but important 
to this section. From now on, we give the following theorem with respect to localization of the QW induced 
by recurrence properties of the corresponding RW. 

Theorem 1 If the RW is positive recurrent with the stationary distribution 7r, then an appropriate choice of 
the initial state gives localization of the corresponding Q W. Moreover, its time averaged limit measure with 
the initial state M/q = X)fc |fc-j|<i Q fe^(fc,j) * s bounded below as follows. 



£°\i,j)>\(aj,* )\ 2 n(iMj). 



(16) 
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Here 



1 



l + C R 

Proof. The corresponding Jacobi matrix is given by 

/ 



So(j) + (l-5o(j)) 



Po ■ --Pj-i 
<7i • • • Qj 



(17) 



y/pTq2 








y/Piq2 







V 



/ 



The detailed balance condition Try^ = Trj+igi+i for z > and positive recurrence of the RW provide J = 
D 1 / 2 MD~ 1 / 2 and Eq. Q7J), where D is a diagonal matrix with (-D)j.j = tt(i). Therefore from Lemma[51 we 
obtain the desired conclusion. 

□ 

Remark 6 If we want to avoid localization of the QW launched at j induced by the positive recurrent RW, 
at least, we should choose the three parameters of initial state atj—i, otj, <Xj+x so that tp = T [aij—i,atj,atj+i] 
is orthogonal to T [y/qj, y/rj, y/Pj\ ■ 

Let us consider a special case rj = for all j € Z + in the following. In this case, from Remark|31 it is obtained 
that % s = 0. Now we check whether —1 is eigenvalue of U or not. In a heuristic argument, a "signed" - 
detailed balance condition Tr'(j)pj = — 7r'(j + 1)^+1 implies that ir'(j) = (— • • -po)/{qj ■ • • qi)- 

Actually, we see 

Mtt' = -iv'. (18) 

Recall that p(x) is the eigenvector of J with its eigenvalue x defined in Eq. ([S]). Combining Eq. (fT5)l with 
Remark 2] gives pj(— 1) = (— l) J ^)j(l) which implies that U has the eigenvalue —1 from Eq. ([7]) in Lemma 
[T] Also we know from Eq. (??) in Lemma [T] that one of the derivations of the value —1 is the eigenvalue of 
H^ s \ but in this case HS S ^ = 0. So we find that the derivation of the eigenvalue —1 comes from just only 
H {R) in this case. Then from Eq. we obtain the following corollary. 



If the RW is positive recurrent, then the time averaged limit measure 

(19) 



Corollary 1 Let rj = for all jsZ 

is bounded below as follows 

Jlg°\i,j)>2x\(a j ,*o}\ 2 n(iMj), 
where ir is the stationary distribution of the RW given in Eq. Ji7[ ). 

Remark 7 Assume that we choose the initial coin state from one of the following set: for j > 1 




with prob. 1/2, 
with prob. 1/2, 



(20) 



for j = 0, \I/o = |0; R). Let us consider the space homogeneous case such that pj — p (j > 1), = 1 (j = 0), 
Qj = Q (i ^ 1)) an d r j — (j > 0) with q > p (positive recurrent). Then the description of section 4 in fllf 
with respect to the presence of two point masses {±1} of the spectral measure of a Chebyshev type stochastic 
matrix, ensures that the time averaged limit measure is described by the lower bound itself in Eq. U9\). that 
is, 

M^ o) (i, j) = {2<5 (j) + (1 - 6 (j)) MjXz), (21) 

where 7r(j) = (1 —p/q)/2 (j = 0), = (1 — p/q)/(2q) x (p/q)^ 1 (j > 1). Thus the summation of the time 
averaged limit measure Cj is given by 



{2*00') 



(1 



S (j))}n(j) < 1, 



which is consistent with the argument of U2\j . which treats a space homogeneous QW on a finite path, in the 
limit of its system size. 



S 



Localization of Theorem Q] and its associated arguments come from H( R ' C %. Recall that the derivation 
of is the corresponding RW. Since the subspace is invariant under the action of U, the choice of 
initial state from T-L^ leads us considerations outside of its orthogonal complement 1-L^ S \ From now on, we 
will face the consideration of %^ s \ For a simplicity, we change the notations of the basis of % as follows: 
S U+hj) -> \j> R )> S (j-i,j) \j> L ), and if r j + °, then -> \ j\0). 



Lemma 3 Assume that S = ■ ■ ■ }, {ji < ji < ■ ••)■ 

- 1 rs: . i>L 

Then it is obtained that T-L^ is spanned by {'r) k }j k £s, where 



+ (-i) Jfc+1 ^ +1 (-^TljWi;i) + \pr~^ l3k+1]0) ) ■ (22) 



iJere = qj (j > 1), =1 (j = 0). In particular, if \S\ — n < oo, i/i 



(7/ 



n n = {-iy-Q Jn (-Jth^h.\j n -,o) + y /q-\j n -,R)) + (-^ l Qi[-Vpi\^ L ) + V9i\hR))- (23) 



In this paper, we call rj fc "signed reflected vector" since f7rj fc = — r/ k . Now we give the proof of Lemma 

m 

proof. Remark that = span{oj ■, bj : j > 0}, where a } = ^fq~j\j;L) + y/rj\j;0) + y/pj\j;R), and 

bj = S aj = ^qj\ 3 -1;R) + ^rj\j; O) + ^pj\j + 1;L). Put e H< s ) as 

= E a 5 r) b'; L > + a 5 0) b'; °> + U; ^ ( 24 ) 

i>o 

where if r*j = 0, then ay' = 0. Then (0, a 3 ) = (<p, bj) = for any j implies 



The above equations lead to 



On the other hand, from the j = cases in Eqs. (|25p . (j2"o| . we have ' = a$ . After all, we get 



a, = a 



= for j > 0. Here we should note that if rj — 0, then — and ay = —-\/q~jJpjOtj-i, 

otherwise, on = —^/qj/rjCtj-i — \pPjJ r ~] 0L j- Substituting these relations into RHS of Eq. (|24[) . we arrive 
at the conclusion. 

□ 
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recurrent 


transient 


|<S| = 





|.S| = 1 





{ji < 3} 


1 <\S\<n 


{h <j< in] 


{ji < 3} 


\S\ = oo 


{h < J} 



Table 2: supp(%( s )) 



For any subspace W C H, define 

supp(H') = { j El+:^<pEU' and 3 J 6 {R, O, L}, s.t. (<£|j; J) ^ 0} , 

where Z+ = {0, 1,2,...}. From Eq. (|T4|) . when we choose the initial state so that \&o G span{q ± (a;o) : 
xq E B}- 1 , then the contribution of localization of the QW is only ||Il5 , 5'o|| 2 , where B is defined in Eq. 
(fl3|) . We should remark that the support of localization is included by supp^H^'), that is supp({Ils , 4'o}) Q 

The subspace T-L^ is outside of Ti.^ which is induced by the corresponding RW. However the recurrence 
properties of the RW appear again roundabout to the following theorem with respect to supp( - H < -' s - ) ). (see 
also Table 2.) 

Theorem 2 Suppose that the set of self loops is denoted by S = {ji, j2, ■ • • }, (ji < 32 < ■ • • ) with jk+i — jk < 
oo. Then it is obtained that 

1. \S\ = case. supp(-H (s) ) = 0. 

2. 1 < \S\ = n < oo case. 



supp("H 



(S)) 



{j E 1i+ : ji < j} if the RW is transient, 

if the RW is recurrent and n = 1, 
{j E Z + : ji < j < j n } if the RW is recurrent and n > 1. 



3. \S\ 



oo case. 



supp(H (s) ) = {j G Z+ :j 1 <j}. 



Proof. If \S\ = 0, then the graph becomes a tree structure. Then from Remark [3l H.^ — which is the 
conclusion of part (1) of the Theorem 2. Next we consider the case |<S| = oo. From LemmalU all the signed 
reflected vectors ||r/ fc || < oo (A; > 0), that is, r] k E £ 2 (A) since supp({?7fc}) = jk+i — jk < oo. Thus we get 
the conclusion of part (3) of Theorem 2. Finally, we consider the case < |<S| = n < oo. In this case, the 
assumption jk+i — jk < oo ensures for all Q < j < n — 1, ||t7j|| 2 < oo. But the £ 2 summability of the final 
signed reflected vector r/ n is not ensured by the assumption because the support of rf n is {j > n}. Now we 
should check the value of ||J7„|| 2 in the following. Equation (|2"3")l in Lemma [3] provides 



E Q 



Note pi + qi = 1 for I > j n , since j n the right most self loop. By using this fact, we get 



(27) 



E Qf= E qi "' qi ~ 1 

2 E 



E 



<7i • --qiPi+qi 



91 ■ ■ ■ 9/ , I /-, r \ 9i " " " 9i» , r 1 - 1i 



-\ oo 



Pi ■ ■ ■ pi 



+ (l-I {j „=o } ) 



Pi •••ft, 



+ / {^=o}- 



Pl 



(28) 



When the RW is transient, then X)tj„+i(9i ' ' ' Qi)/(Pi "'Pi) < Ct < oo. Combining Eqs. (|2"T|) with 
we obtain the desired conclusion of part (2). Then we complete the proof. 
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□ 

In the following two corollaries, to highlight effects of to localization, we impose the following assump- 
tions to the initial state \IV 

(1) * € span{q ± (x ) : x G B}^, 

(2) supp(* ) < oo. 

Define ~px°\j) = h m T->oo ^/TYin=n P(X^ ^ = j). We give simple examples of Theorem 2 as the following 



corollaries which give time averaged limit measures by directly computing 



J£{L,0,R} 



l(*;J|n 5 *c 



Corollary 2 Suppose that S = {0}. If the RW is transient, then the time averaged limit measure of the 
corresponding QW with the initial state ty a satisfying the assumptions (1) and (2) is given by 



where 



a , 



j>0 



1 r^W + d-W)^^-- 1 
1 Pi---Pj 



l + O 



R 



V ^ V V; l-^|0,O)+v^|0,i2) :j=0. 



(29) 

(30) 
(31) 



Here 



In the case of a RW, the existence of finite number of self loops have little influence to the behavior of the 
RW itself. However as is seen in the next corollary, the existence of only one self loop changes localization 
property of the corresponding QW. 

Corollary 3 In a recurrent RW with S = {0}, localization of the corresponding QW does not occur for any 
initial states in 'J/o satisfying the assumptions (1) and (2). On the other hand, once we add an additional 
self loop at position n > 0, then an appropriate choice of initial state from span{q ± (a;o) : xo G B}- 1 so that 
|(a^'^ n \*o)| > gives localization whose time averaged limit measure is described by 



where 



^<«j- (B \*o> 



j>0 



7T' (n) «, 



(1 

ho<3<n} x I r 



roqi--- qj -i j 

Pl - Pj X J {">1} 

I mi^zi x (i- Pn ) 



J = o, 

l<j<n-l, 

3 = n, 



, v^|0,O) + V^|0,i?) 

4' (n) = (-iy^n' in) (j) x { {-JPj\j,L) + Jqj\j,R)) x I {n>1} 

y/pH\n, L) + ^Jp n q n /r n \n, O) 



j = 0, 

1 < j < n - 1, 
j = n. 



(32) 



(33) 



(34) 



Here 



C { n y = I {n >i } x{J2 



roqi---qj-i r qi---q n -i x n N [ , T r (l-pi) 



+ 



—[ Pi---Pj Pi---Pn-ir n 



X (1 -Pn) ) +I{n=l} : 



n 



n 



We have found that the time averaged limit measure is described by the summation over all the orthogonal 
projections of the initial state onto each £ 2 (A) eigenvector of U. The £ 2 (V) eigenvectors of the corresponding 
RW give all the £ 2 (A) eigenvectors in In general, it is difficult to obtain all the £ 2 (V) eigenvectors 

except the maximal eigenvector corresponding to the stationary measure of the RW. By using the available 
eigenvectors (including the stationary measure) , we have given a lower bound of the time averaged limit mea- 
sure and evaluate localization of the QW. In some cases, all the £ 2 (V) eigenvectors of the RW corresponding 
to all the mass points of the spectral measure can be computed. For example, from a quantum probabilistic 
approach, in the case of pj = p, qj = q and rj = r with pj + qj + rj = 1, the spectral measure of the Jacobi 
matrix is expressed by the free Meixner law [13J . In Ref. |14j . by using this, the time averaged limit measure 
in this case is explicitly obtained. Moreover they give an expression for a complete orthonormal system of 
in a structure of the Sierpinski gasket, and then a condition of the initial state for localization with a 
bounded support can be shown. 

4 Discussion 

We decomposed the eigensystem of the QW into induced by the RW and its orthogonal complement 

%( s > whose eigenvalue is —1. We showed that the fundamental derivation of localization is the orthogonal 
projection of the initial state onto I 2 {A) eigenvectors which derives from £ 2 (V) eigenvectors of the RW, and 
also onto the subspace W s \ It was obtained that sufficient conditions for localization are expressed by 
recurrence properties of the RW. These are quenched results. 

Due to the infiniteness of the system, the distributivity of the time averaged limit measure is not ensured. 
Let Co be the summation of the time averaged limit measure over all the positions. Then how can we 
characterize the value 1 — c$? From now on, we call the value 1 — cq "missing value". For example, when 
Pj = P (j > 1)) = 1 (j = 0), and qj — q (j > 1) with p + q = 1, then according to [15], we see that the 
missing value appears at the following weak limit theorem. 

^±^c 5 (x) + (l-c )- X - rf K (x:2^pq) (t -> oo), (35) 

t l-\p-q\ 

where Jk{x : s) is the Konno density function [TH [T7] with parameter s € (0, 1), cq — I{ P < q }(^ — p/q) and 
means the weak convergence. This is so called "clean" quantum walk case. Now how about disordered 
quantum walk case? In the case of one defect on the doubly infinite line [18], the missing value, which 
spreads away ballistically, appears in the weak limit theorem which is expressed by a similar formulas of Eq. 
(|35j) : the density is a convex combination of a delta measure at the origin corresponding to localization, and 
a product of a rational polynomial and the Konno density function. This result of [T5] suggests that a small 
number of perturbations generates edge state of the QW, and removes the missing value corresponding to 
bulk state ballistically to infinite places keeping a shape characterized by the Konno density functions as 
is the second term of RHS of Eq. ([35]) . On the other hand, interestingly, Obuse and Kawakami [19] and 
Ahlbrecht et al. [20j suggest that a large number of perturbations, that is, an random environment, gives 
a characterization of the missing value by an exponential function named dynamical localization with not 
ballistic but sub-ballistic spreading. 

Finally, the quantum coins treated in this paper belong to a special class of unitary operator. The 
stochastic process behind this QW with this choice of quantum coins makes the computations be easy and 
gives a nice relationship of behaviors between the RW and the QW. What kind of walk exists behind the 
QW or the CMV matrix [21] with another choice of quantum coins such as [22]? Applying the frame work 
in this paper to another kind of class of quantum coins is one of the interesting future problems. 
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